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Abstract 

The nonleading corrections in the inverse mass of the heavy charmed quark to 



a* 



X 



the effective local AS* = 2 Lagrangian of K° — K° mixing are calculated within the 
standard model. 
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1. Nowadays the standard model of strong and electroweak interactions [1-3] is 
the subject of intensive and thorough experimental study [Q. So far the model goes 
successfully through the numerous precious tests even though the accuracy of the 
experimental data is such that the radiative corrections must be properly taken into 
account [[|. One of the most interesting effects having got the explanation within the 
standard model is the phenomenon of CP-invariance breaking observed firstly in the 
system of the neutral kaon but later found for the B° — B° as well. The problem of 
CP-invariance breaking is still very intriguing and both experimental and theoretical 
investigations of the K° — K° system draw much attention. In this particular case, 
however, the direct comparison of the theoretical predictions with experimental data is 
made somewhat difficult due to the necessity of determination of the hadronic matrix 
elements of the effective electroweak Lagrangian. It is this problem which causes the 
main part of uncertainties of the theoretical calculations. The effective electroweak 
local AS = 2 Lagrangian for the K° — K° mixing in the leading approximation in the 
heavy charmed quark mass m c has been obtained in the famous paper by Gaillard and 
Lee [[J. Then the leading perturbation theory corrections due to strong interactions 
were calculated in 0. The corresponding hadronic matrix element of the effective 
local L&s=2 Lagrangian between K° and K° states 

(K°(k')\L AS=2 \K°(k)) 

has been intensively studied during several last years with different techniques [8-17] 
and, probably, the final result is consistent with the vacuum saturation hypothesis 
within 20% accuracy. 

However, Wolfenstein |1| has pointed out that the local effective Hamiltonian 



does not exhaust the whole physics of AS = 2 transitions. It cannot account for 
the long distance contribution which is present in the initial Green's function for the 
matrix element of the K° — K° mixing and is connected with the propagation of 
the light u— quark round the loop of the box diagram. This contribution is purely 
non-perturbative and ultimately depends on the infrared properties of QCD. There 
are some model dependent estimates of this part of the entire matrix element in the 
literature [19-23] though the accuracy is still far from being satisfactory. 

Nevertheless, there is one point which can be essentially improved just within 
perturbation theory for the standard model. It consists in the calculation of the 
corrections in the inverse mass of the charm quark to the local part of the effective 
Lagrangian. These correction can be important because the charmed quark is not 
heavy enough in comparison with the characteristic scale in the sector of light (u, d) 
quarks, for example, with the p— meson mass, m p . The analogous corrections to the 
decays of charmed baryons and mesons have been recently considered in refs. [24-27]. 

In the present paper the non-leading corrections in m" 1 to the local part of the 
effective AS* = 2 Lagrangian are calculated. They are represented by the local oper- 
ators with dimension eight in mass units. 

2. The local effective AS = 1 Lagrangian after decoupling of the W-boson has 



the form 

las=i = ^j,j; (i) 

where Gp is the Fermi constant, J^ = QhlpVqL is the weak charged hadronic current, 
Q = (u, c, t) T , q = (d, s, b) T , q^ denotes the left handed Dirac spinor, and V is the 
quark mixing matrix. Then the matrix element M of the K° — K° transition is 
represented by the expression 

out (K°(k')\K°(k)) m = i(2n) 4 5(k - k')M, 

M = l -J(K°(k')\TL AS=1 (x)L AS=1 (0)\K%k))dx. (2) 

Strictly speaking, formulae ([J) and @ are literally valid for the situation where 
the t— quark is much lighter than the W-boson. As is well known now this is not the 
case anymore. For our purpose, however, it is not essential and in the following we 
neglect the t— quark mixing at all and restrict ourselves to the simplified model with 
two generations only. 

The effective AS = 2 Lagrangian can be written in the form 

L^=( ^^f^ j(L H + L L ) (3) 

where 6 C is the Cabibbo angle and we have used the obvious normalization. Here 

L H = i J T H (x)dx, T H = T cc - T cu - T uc (4) 

is the heavy part of the whole effective AS = 2 Lagrangian containing the loops with 
a virtual heavy c— quark in the intermediate state, while 

L L = i T L (x)dx, T L = T uu (5) 

describes the light part of the transition. We also introduce the useful notations 

T uu (x) = Ts L '~<( a u L u L '~<( a d L (x)s L '~i l 3U L u L '~<ipd L (®), 

T cu (x) = Ts L ^ a u L c L ^ a d L (x)s L ^pc L u L ^f3d L (Q) 

and so on. Expression (3) is finite due to GIM mechanism but Lh and L^ separately 
require a regularization because they are ultraviolet divergent. The dimensional reg- 
ularization is not convenient in this case due to the presence of the 75-matrix. The 
Pauli-Willars regularization introduces a regulator mass that makes it difficult to 
perform an explicit calculation from the technical point of view. We will use the reg- 
ularization which is free of the shortcomings of both of them in our particular case. 
Namely, let us define the regularized quantities L^ L by the equation 

Lh, l = tjT H , L {x){-ii 2 x 2 ydx (6) 



where e is the regularization parameter and \x represents the mass scale analogous to 
one of the dimensional regularization. 

Now for the heavy part of the effective AS = 2 Lagrangian Lh we have a regular 
expansion in the inverse charmed quark mass in the following form (from now on we 
omit the index "R") 

lQn 2 L H = C Oo + J2C 3 Oj 
j 

where Cj are coefficient functions depending on the heavy quark mass m c and Oj are 
the local operators built from the light (u, d, s) quark field only. If we split the whole 
Lagrangian into the sum 

Lh = L<$ + lU (7) 

then 

16n 2 LP = -m 2 c (s Lla d L ) 2 (8) 

is the well known result of Gaillard and Lee . The rest part of eq. (0) 

16^^ = ^0,0, (9) 

i 

consists of local operators which have dimension eight in mass units. It is convenient 
to define the operator basis in the form 

Op = SLjadLSLJuFpadL, 

Oa = s L ^^D v )d L s L ^^D v )d Ll 
B = s L 'j ll D f ,d L s L 'y^D^d L , 

( 771 — \- 777 I 
Oc = S L "/ a d L SL( r rn D » D a + DoT/^D^dL S — —(s L lad L ) 2 

where 



IfQiDu) = -fopDv+lvDp)- 



Then a direct calculation gives 



^-^ + o A) I + **£:,+§ -| 



- ^{0 B + O c ) ( - + In (-^—) + - ) (10) 



2 ^ r, ^( 1 , Me- 20 , 11' 
77 (0 B + O c - + hi + — 

3 ye m 2 c 6 / 

where C = 0.577 ... is the Euler constant. 

Eq. ( |10"D is the main result of the present paper. After performing a renormaliza- 
tion procedure (say, a minimal subtraction of the pole term) we will have the finite 
quantity and the parameter \x recalls the necessity to have the proper short distance 
contribution missed in eq. (0) and given by eq. (f|). At this order of the expansion 



in m" 1 the dependence on this parameter is explicit, contrary to the leading order 
(|8|) which is finite and does not depend on /x at all. It is obvious from the expres- 
sion (|T0|) that the heavy and light parts of the whole Lagrangian must be defined 
simultaneously in a coordinated way. The pole part of eq. fllDD is canceled by the cor- 
responding divergences of the light part due to GIM mechanism. The last statement 
can be explicitly demonstrated. 

Namely, let us consider the light part. The operator product expansion for the 
amplitude T L (x) from eq. fl5|) in x 2 at x 2 — > has the form 

T L {x) = ^(s Lla d L ) 2 + Ij ^- J (±{pz + A ) + \{0 B + O c )) ■ (11) 

The short distance contribution of the light part does cancel the corresponding 
divergences of the heavy part. More technically we extract the short distance con- 
tribution of the light part as follows. We split the entire light part in the following 

way 

Ll = i TL(x)(—fi 2 x 2 ) £ dx 

= ijT L {x){f{x,x) + f(x,x))(-ii 2 x 2 ) e dx = L S L H + L L L G (12) 

where 

L S L H = i JT L {x)f{x,x){-n 2 x 2 fdx, (13) 



and 



L LG = % \ T L (x)f(x, x)(-fi 2 x 2 Ydx, (14) 



f(x,x) + f(x,x) 



A simple model which is the smooth generalization of the step #— function is chosen 
for the functions f(x, x) and f(x, x) 



rv* £> Th ry ■" ■" 



f^ S) = ^ , ? ,,„ , f(x, X 



x 2n + (-x 2 ) n ' X 2n + (-x 2 ) n 

in such a way that the function f(x,x) cuts out the short distances only (up to x) 
and the function f(x, x) does the long distances. After performing the integration we 
get for the short distance contribution of the light part the expression 

167r 2 Lf = (s LlQ d L ) 2 -^-- 2 + (\(Oz + O a ) + \{0 B + O c j) (- + \n^ 2 x 2 

sm(7r/nj x l V3 3 / \e 

Now the whole answer is 



16tt 2 (L h + L L ) = 1Qtt 2 L^ + (s Lla d L ) 2 \-m 2 c + - 



sin(7r/n) x 



7-2 



^(o P + o A ) + ^(o B + o c ; 



'"W + 3, 



2 -o A . 

3 



(15) 



The long distance part of the whole light Lagrangian L£ ([14] ) cannot be calculated 
due to strong infrared problems and requires some low-energy model. It can be 
estimated from lattice calculations or with the help of the theory of effective chiral 
Lagrangians. We hope to consider this question in a separate publication. 

To get the numerical estimates for the corrections obtained in the paper (eq. (|i~5])) 
we use factorization for the kaon-antikaon matrix elements of the local operators 



Op, — Oc- In the chiral limit we have 



(K°(k)\(s L 7ad L ) 2 \K°(k)) 



fact 



^J 2 K m 2 K . 



N, 



)( : 



{K (k)\O F \K°(k)) 



fact 



.j%m 2 

2 



A' 



{K (k')\O A \K°(k)) fact = S 2 ^-(^^ 



N r 



where the parameter 5 is defined by the relation 



\sLl^ a d L \K°{k)) = -ik^f K 5 2 



and all other matrix elements vanish. 
From eq. ([15]) we finally have 



16tt 2 (L h + L b L 



SH\ 



(K\k)\{s Lla d L ) 2 \K\k)Y act \-m. 
4, 1 



n/n 



sin(7r/n) x 



7-2 



4 r2„ ,*e- 2C 



+ 



6 



-5 2 



Let us estimate x 2 . From eq. (|TT|) we find 
1 



2^2 



S 2 x 



I, (,-) = ^^(K°(k)\(s Lla d L ) 2 \K°(k)y^(l - u -f- + o(x 2 )). 



(16) 



(17) 



Then at 5 2 x 2 = 3 where the expansion ( |i~7|) blows up and at n = oo we obtain the 
following representation for the non-leading corrections 

-ml < K\k)\{s L l a d L ) 2 \K\k) >/- (l-l^-l^(ln(^)-2C + |)). (18) 



Numerically at m 2 c = 1.6 GeV 2 , 6 2 = 0.2 GeV 2 [||, || we get 
-m 2 c (K°(k)\(s Lla d L ) 2 \K°(k)y act (l - 0.2 + 0.3) 



= -m 2 c {K\k)\(s Lla d L f\K\k))f ac \l + 0.1). 

But if we define the boundary value of x 2 by the relation S 2 x 2 = 1 then the 
corrections are much larger and have the opposite sign 

-m 2 c (K°(k)\(s Lla d L ) 2 \K°(k)y act (l - 0.5 + 0.1) 

-m 2 c (K°(k)\(s Lla d L ) 2 \K°(k)) fact (l - 0.4). 

3. To conclude, we have calculated the corrections in the inverse mass of the 
charmed quark, considered as heavy enough, to the effective local AS" = 2 Lagrangian. 
In this order there already exists the explicit dependence on the boarder between short 
and long distances even in the heavy part of the whole Lagrangian. A splitting of 
light part into the short and long distance contributions is proposed which respects 
the cancelation of purely ultraviolet divergences, i.e. the GIM mechanism. The OPE 
at x 2 — > for the light part of Lagrangian is analyzed and the convergence scale is 
determined explicitly. In the vacuum dominance approximation this scale is large 
enough and is determined by the parameter 5~ 2 . The corrections can decrease the 
answer for the transition matrix element approximately twice though it depends on 
the concrete choice of the scale x 2 and also on the validity of the vacuum dominance 
approximation for estimating the kaon matrix elements of the operators F — Oq- 
This means that the purely long distance hadronic contribution could be important to 
stabilize the whole answer in the region where the expansion (13) fails, in other words 
one can predict that the integral Ljf* changes quickly when x 2 belongs to the region 
1/5 2 < x 2 < 3/5 2 just to compensate the corresponding dependence of the short 
distance part flilf ). Nevertheless, in ordert to determine the absolute magnitude of 
the long distance hadronic contribution one needs a low-energy model such as, for 
example, lattice or chiral Lagrangians. 
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